Vortex formation in a stirred Bose-Einstein condensate 
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Superieure 



Using a focused laser beam we stir a Bose-Einstein con- 
densate of 87 Rb confined in a magnetic trap and observe the 
formation of a vortex for a stirring frequency exceeding a crit- 
ical value. At larger rotation frequencies we produce states 
of the condensate for which up to four vortices are simulta- 
neously present. We have also measured the lifetime of the 
single vortex state after turning off the stirring laser beam. 
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Rotations in quantum physics constitute a source of 
counterintuitive predictions and results as illustrated by 
the famous "rotating bucket" experiment with liquid he- 
lium. When an ordinary fluid is placed in a rotating con- 
tainer, the steady state corresponds to a rotation of the 
fluid as a whole together with the vessel. Superfluidity, 
first observed in liquid Hell, changes dramatically this 
behavior [QJ|]. For a small enough rotation frequency, 
no motion of the superfluid is observed; while above a 
critical frequency, lines of singularity appear in its ve- 
locity field. These singularities, referred to as vortex 
filaments, correspond to a quantized circulation of the 
velocity (nh/m where n is an integer, and m the mass 
of a particle of the fluid) along a closed contour around 
the vortex. In this letter we report the observation of 
such vortices in a stirred gaseous condensate of atomic 
rubidium. We determine the critical frequency for their 
formation, and we analyze their metastability when the 
rotation of the confining "container" is stopped. 

The interest in vortices for gaseous condensates is that, 
due to the very low density, the theory is tractable in 
these systems and the diameter of the vortex core, which 
is on the order of the healing length, is typically three 
orders of magnitude larger than in Hell. At this scale, 
further improved by a ballistic expansion, the vortex fil- 
ament is large enough to be observed optically. The 
generation of quantized vortices in gaseous samples has 
been the subject of numerous theoretical studies since 
the first observations of Bose-Einstein condensation in 
atomic gases ||-||. Two schemes have been considered. 
The first one uses laser beams to engineer the phase of 
the condensate wave function and produce the desired 
velocity field [0-0- Recently this scheme [[U)| has been 
successfully applied to a binary mixture of condensates, 
resulting in a quantized rotation of one of the two com- 



ponents around the second one |12|. Phase imprinting 
has also been used for the generation of solitons inside a 
condensate 

The second scheme, which is explored in the present 
work, is directly analogous to the rotating bucket ex- 
periment |T^ , |l6| ]. The atoms are confined in a static, 
cylindrically-symmetric Ioffe-Pritchard magnetic trap 
upon which we superimpose a non-axisymmetric, attrac- 
tive dipole potential created by a stirring laser beam. The 
combined potential leads to a cigar-shaped harmonic trap 
with a slightly anisotropic transverse profile. The trans- 
verse anisotropy is rotated at angular frequency fl as the 
gas is evaporatively cooled to Bose-Einstein condensa- 
tion, and it plays the role of the bucket wall roughness. 

In this scheme, the formation of vortices is -in 
principle- a consequence of thermal equilibrium. In the 
frame rotating at the same frequency as the anisotropy, 
the Hamiltonian is time-independent and one can use 
a standard thermodynamics approach to determine the 
steady-state of the system. In this frame, the Hamil- 
tonian can be written H = H — £IL Z , where H is the 
Hamiltonian in the absence of rotation, and L z is the 
total orbital angular momentum along the rotation axis. 
Above a critical rotation frequency, f2 c , the term —VLL Z 
can favor the creation of a state where the condensate 
wave function has an angular momentum h along the 
z axis and therefore contains a vortex filament [p7|-p5[. 
The density of the condensate at the center of the vortex 
is zero, and the radius of the vortex core is of the order 
of the healing length £ = (87rap) -1 / 2 , where a is the scat- 
tering length characterizing the 2-body interaction, and 
p the density of the condensate |2(| . 

The study of a vortex generated by this second route 
allows for the investigation of several debated questions 
such as the fate of the system when the rotating velocity 
increases above f2 c . This could in principle lead to the 
formation of a single vortex with n > 1 at the center of 
the trap; however, this state has been shown to be cither 
dynamically or thermodynamically unstable ]l|pl|- p^ , p7t . 
The predicted alternative for large rotation frequencies 
consists of a lattice of n — 1 vortices. Another important 
issue is the stability of the current associated with the 
vortex once the rotating anisotropy is removed p8| . 

Our experimental set-up has been described in de- 
tail previously [^9|. We start with 10 9 87 Rb atoms in 
a magneto-optical trap which are precooled and then 
transferred into an Ioffe-Pritchard magnetic trap. The 
evaporation radio frequency starts at u T { — 15 MHz 
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and decreases exponentially to i/^f na1 ' in 25 s with a 
time constant of 5.9 s. Condensation occurs at Av^ = 
^(finai) _ ^ 2 g ^O 6 atoms and a tem- 

perature 500 nK. Here = 430 (±1) kHz is the radio 

frequency which empties completely the trap. The slow 
oscillation frequency of the elongated magnetic trap is 
lu z /(2tt) — 11.7 Hz (z is horizontal in our setup), while 
the transverse oscillation frequency is lj±/(2tt) = 219 Hz. 
For a quasi-pure condensate with 10 atoms, using the 
Thomas-Fermi approximation, we find for the radial and 
longitudinal sizes of the condensate Aj_ = 2.6 /im and 
A z = 49 /im, respectively. 

When the evaporation radio frequency reaches the 
value i/ r ( f min) + 80 kHz, we switch on the stirring laser 
beam which propagates along the slow axis of the mag- 
netic trap. The beam waist is w s = 20.0 (± 1) /tm and 
the laser power P is 0.4 mW. The recoil heating induced 
by this far-detuned beam (wavelength 852 nm) is negli- 
gible. Two crossed acousto-optic modulators, combined 
with a proper imaging system, then allow for an arbi- 
trary translation of the laser beam axis with respect to 
the symmetry axis of the condensate. 

The motion of the stirring beam consists in the super- 
position of a fast and a slow component. The optical 
spoon's axis is toggled at a high frequency (100 kHz) 
between two symmetric positions about the trap axis z. 
The intersections of the stirring beam axis and the z = 
plane are ±a(cos9u x + sin6*Uj,), where the distance a is 
8 /im. The fast toggle frequency is chosen to be much 
larger than the magnetic trap frequencies so that the 
atoms experience an effective two-beam, time averaged 
potential. The slow component of the motion is a uni- 
form rotation of the angle 8 — Sit. The value of the 
angular frequency SI is maintained fixed during the evap- 
oration at a value chosen between and 250 rad s _1 . 

Since w s 3> Aj_, the dipole potential, proportional to 
the power of the stirring beam, is well approximated 
by mLU 2 ± (e x X 2 + e Y Y 2 )/2. The X, Y basis is rotated 
with respect to the fixed axes {x,y) by the angle 0(t), 
and ex = 0.03 and ey = 0.09 for the parameters given 
above |50j . The action of this beam is essentially a slight 
modification of the transverse frequencies of the mag- 
netic trap while the longitudinal frequency is nearly un- 
changed. The overall stability of the stirring beam on the 
condensate appears to be a crucial element for the suc- 
cess of the experiment, and we estimate that our stirring 
beam axis is fixed to and stable on the condensate axis 
to within 2 /im. We checked that for SI < Sl c the stirring 
beam does not affect the evaporation. 

For the data presented here, the final frequency of the 
evaporation ramp was chosen just above i/f" (Ai^f G 
[3,6] kHz). After the end of the evaporation ramp, we 
let the system reach thermal equilibrium in this "rotating 
bucket" for a duration t r = 500 ms in the presence of an rf 
shield 30 kHz above uS nal ^ . The vortices induced in the 



condensate by the optical spoon are then studied using a 
time-of- flight analysis. We ramp down the stirring beam 
slowly (in 8 ms) to avoid inducing additional excitations 
in the condensate, and we then switch off the magnetic 
field and allow the droplet to fall for r = 27 ms. Due to 
the atomic mean field energy, the initial cigar shape of 
the atomic cloud transforms into a pancake shape during 
the free fall. The transverse xy and z sizes grow by a 
factor of 40 and 1.2 respectively [[h). In addition, the 
core size of the vortex should expand at least as fast as 
the transverse size of the condensate 
a vortex with an initial diameter 2£ 
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-33]. Therefore 
0.4 /jm for our 
experimental parameters is expected to grow to a size of 
16 /im. 

At the end of the time-of-flight period, we illuminate 
the atomic sample with a resonant probe laser for 20 /is. 
The shadow of the atomic cloud in the probe beam is 
imaged onto a CCD camera with an optical resolution 
~ 7 /im. The probe laser propagates along the z-axis so 
that the image reveals the column density of the cloud 
after expansion along the stirring axis. The analysis of 
the images, which proceeds along the same lines as in 
p9| , gives access to the number of condensed Nq and un- 
condensed N' atoms and to the temperature T. Actually, 
for the present data, the uncondensed part of the atomic 
cloud is nearly undetectable, and we can only give an 
upper bound for the temperature T < 80 nK. 

Figure [j] shows a series of five pictures taken at vari- 
ous rotation frequencies SI. They clearly show that for 
fast enough rotation frequencies, we can generate one or 
several (up to 4) "holes" in the transverse density dis- 
tribution corresponding to vortices. We show for the 0- 
and 1-vortex cases a cross-section of the column density 
of the cloud along a transverse axis. The 1-vortex state 
exhibits a spectacular dip at the center (up to 50 % of 
the maximal column density) which constitutes an un- 
ambiguous signature of the presence of a vortex filament. 
The diameter of the vortex core following the expansion 
is measured at the half max of the dip to be ~ 20 /im. 

For a systematic study of the vortex stability domain, 
we have varied in steps of 1 Hz the rotation frequency 
for a given atom number and temperature. For each fre- 
quency, we infer from the absorption image the number 
of vortices present, and the results are shown in Fig. Q. 
Below a certain frequency, we always obtain a conden- 
sate with no vortices. Then, in a zone with a 2 Hz width, 
we obtain condensates showing randomly or 1 vortex. 
Increasing SI, we arrive at a relatively large frequency in- 
terval (width 10 Hz) where we systematically observe a 
condensate with a single vortex present. Our value for 
the critical frequency is notably larger than the predicted 
value of 91 Hz |§ (see also [0jllj2|,|2| ) . This devia- 
tion may be due to the marginality of the Thomas-Fermi 
approximation for our relatively low condensate number. 
If SI is increased past the upper edge of the 1-vortex zone, 
multiple vortices, as shown in Fig. |l], are observed . The 
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range of stability of the multiple vortex zones appears to 
be much smaller than that for the 1-vortex zone. The 
3-vortex zone, for instance, seems to be stable over only 
3-4 Hz and is complicated by the occasional appearance 
of a 2-vortex or a 4-vortex condensate. At this stage of 
the experiment, it is difficult to determine whether these 
shot-to-shot fluctuations are due to a lack of experimental 
reproducibility or to the fact that these various states all 
have comparable energies and therefore all have a reason- 
able probability to occur for the range of parameters in 
question. Finally, when f2 is increased past the range of 
stability for the multiple vortex configuration, the density 
profile of the condensate takes on a turbulent structure, 
and the condensate completely disappears for f2 larger 
than 210 Hz, which should be compared with the aver- 
age transverse frequency of the magnetic + laser dipolc 
potential (226 Hz). 

It is remarkable that the multiple vortex configurations 
most often occur in a symmetric arrangement of the vor- 
tex cores: an equilateral triangle and a square for the 
3-vortex and the 4-vortex cases respectively. This find- 
ing supports the theoretical analysis of (24|] which shows 
that vortices rotating in the same direction experience an 
effective repulsive interaction, which in turn favors these 
stable configurations (see also |pl"| ). 

The final question addressed in this letter concerns 
the lifetime of a vortex state in an axisyrrrmetric trap. 
Without a rotating anisotropy, the vortex state is no 
longer the lowest energy state of the system, and after 
the anisotropy is removed, one expects that the gas will 
eventually relax to a condensate with no vortex plus a 
slightly larger thermal component, bolstered by the en- 
ergy contained in the vortex state. Figure [| presents the 
experimental study of the single-vortex state lifetime at 
two different condensate parameters. We choose a rota- 
tion frequency f2 in the middle of the 1-vortex range of 
stability, and we let the vortex form as before in the pres- 
ence of the stirring beam. Then we switch off the stirring 
beam, and we allow the gas to evolve in the pure mag- 
netic trap for an adjustable time. Finally, we perform 
a time-of-flight analysis to determine whether the vortex 
is present or not. Each point in the two lifetime curves 
represents the average of 10 shots, where we have plotted 
the fraction of pictures showing unambiguously a vortex 
as a function of time J34|. We deduce from this curve a 
characteristic lifetime of the vortex state in the range 400 
to 1000 ms with a clear non-exponential decay behavior. 
In addition, we observed that at long times the vortex 
rarely appears well centered as it does immediately after 
formation. 

To summarize, we have reported the formation of vor- 
tices in a gaseous Bose-Einstein condensate when it is 
stirred by a laser beam which produces a slight rotating 
anisotropy. A natural extension of this work is to study 
the superfluid aspects of this system j3?| by investigating 
the dynamics for vortex nucleation and decay. An impor- 



tant question is the role of the thermal component. For 
instance, nucleation can occur either by transfer of an- 
gular momentum from this component to the condensate 
or directly from a dynamical instability of the non- vortex 
state p2] , p4] , |3"6| , |37| . Also the decay of the single vortex 
state may be due to the coupling with a non-rotating 
thermal component fl38| or due to the instability induced 
by the residual, fixed anisotropy of our magnetic trap 
(measured to be u x /uj y = 1.012±.002) ||§. Finally, this 
type of experiment gives access, in principle, to the el- 
ementary excitations of the vortex filament pf]|-^2[|, the 
study of which might reveal new aspects of the superfluid 
properties of these systems. 
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FIG. 1. Transverse absorption images of a Bose-Einstein 
condensate stirred with a laser beam (after a 27 ms 
time-of-fiight) . For all five images, the condensate number 
is No = (1.4 ± 0.5) 10 5 and the temperature is below 80 nK. 
The rotation frequency 0/(2-71") is respectively (c) 145 Hz; (d) 
152 Hz; (e) 169 Hz; (f) 163 Hz; (g) 168 Hz. In (a) and (b) we 
plot the variation of the optical thickness of the cloud along 
the horizontal transverse axis for the images (c) (0 vortex) 
and (d) (1 vortex). 



FIG. 2. Edges of the regions of stability for the 0, 1, and 
multiple vortex configurations. The condensate number is 
N = (2.3 ± 0.6) 10 5 and the temperature below 80 nK. 



FIG. 3. Fraction of images showing a vortex as a func- 
tion of the time spent by the gas in the axisymmetric trap 
after the end of the stirring phase for two condensate condi- 
tions. The sets of data correspond to condensate numbers of 



N, 



(°) 



= (2.3 ± 0.6) 10 b and A^' 5 = (1.2 ± 0.3) 10 b which 
were obtained with Af,f =6 and 3 kHz respectively. This 
implies that > T^*\ where both temperatures are below 
our detection limit of 80 nK. 
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